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Abstract 

The SO (3) lattice gauge theory with a Villain form of action was 
investigated by Monte Carlo techniques on asymmetric lattices with 
iV T = 2 and 4, where N T is the number of sites in the temporal ex- 
tent. Unlike the results for higher iV T , only one transition of second 
order was found for N T = 2 . An extended action with an irrelevant 
term to suppress Z2 monopoles enabled us to get a better view of the 
deconfinement transition as the effects of bulk transition could be sup- 
pressed as well. Although the action has no global ^-symmetry for 
the SO (3) theory, unlike the SU(2) theory at finite temperature, our 
study revealed a second order deconfinement transition, with proper- 
ties similar to the deconfinement transition of SU(2). 
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1 INTRODUCTION 



The most convenient way of regularizing gauge theories for nonperturbative 
studies is by discretizing the underlying space-time on a lattice. The gauge 
variables can then be chosen as suitable representations of the gauge group 
residing on the links of the lattice. 

A natural choice for the lattice action is the Wilson action [I], in which 
the gauge variables are in the fundamental representation of the gauge group. 
For the SU(2) gauge theory, the action is simply 

z v 

where U p is the product of link variables around a plaquette, and the sub- 
script / denotes the trace is to be taken in the fundamental representation of 
the gauge group. The fundamental representation is particularly convenient 
when one includes fermions, since the fermions are normally taken in this 
representation. The Wilson action has been widely studied in the literature, 
both for zero-temperature and finite temperature studies. Studies of gauge 
theories at finite temperature can be performed by taking a lattice with a fi- 
nite temporal extent equal to the inverse temperature, and imposing periodic 
boundary conditions on the gauge variables in the temporal direction. One 
interesting prediction that has emerged from such studies is that, at a suffi- 
ciently high temperature, the theory undergoes a deconfinement transition. 
The finite temperature gauge theory has a global symmetry corresponding 
to the center of the gauge group. The deconfinement transition is gener- 
ally described by the breaking of this center symmetry. The Polyakov loop 
operator, 

Lf(f)=Trf l[U T (f,i), (2) 

%=i 

transforms non-trivially under this symmetry and therefore acts as an order 
parameter for the deconfinement transition. Here N T is the number of lattice 
points in the temporal direction. 

The Wilson action (P is by no means the only choice for the discretized 
action. Since the continuum theory is at a critical point of the lattice regu- 
larized theory, one expects a large class of actions, differing only by irrelevant 
terms, to correspond to the same continuum physics. In particular, one can 
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take the trace in (Q) in any representation of the gauge group. Since the con- 
tinuum theory is defined only in terms of the underlying algebra, the choice 
of the representation for the link variables should be completely irrelevant 
for it. 

One particularly interesting variant of ([[]) is when the trace is taken in 
the adjoint representation. The gauge variables in adjoint representation are 
blind to the center of the group. The center symmetry is therefore absent 
and (Lf) is identically zero. A study of the deconfinement transition for this 
form of the action is therefore of particular interest ||. 

Finite temperature 5*0(3) gauge theory was studied in refs. || f| with 
the Wilson form for the action, and in ref. || using a Villain form ||. A 
natural analog of (Lf) was taken to be (L a ) where the adjoint Polyakov loop 
L a is defined as 



A'r 



L (f)=Tr Q I[U r {r,i). (3) 

While (L a ) is not an order parameter for the deconfinement transition as the 
screening of an adjoint quark renders it to be non-zero in the confined phase, 
it is known to behave quite similarly as (Lf) across the SU(2) deconfinement 
transition. It was found in refs.[]3], [|, [5| that the finite temperature theory 
shows a discontinuous transition. (L a ) ~ till the transition point, where it 
rises sharply, indicating a deconfining nature of the transition. The transition 
point, however, was coincident with the known bulk transition point for 
5*0(3) 0. Moreover, it did not shift with increasing temporal extent of the 
lattice [|], |5j, which is both unlike the deconfinement transition in the SU(2) 
theory and characteristic of a bulk transition. 

These results appear to be similar to those obtained in refs. || |9| in finite 
temperature studies with the Bhanot-Creutz action |F) 



s = £ [P/ (i - ^ Tr /^) + & (i - l T ^ u v) 



(4) 



There too, a first order deconfinement line, as characterized by the order 
parameter (Lf) , was found to be coincident with the previously known first 
order bulk transition line for a range of large (3 a . Moreover, the shift in 
the transition coupling with changing N T was found to be negligibly small 
although it increased progressively as (3 a was decreased. 

Both for the intrinsic importance of the study of the deconfinement transi- 
tion for the 50(3) gauge theory, and in order to have a better understanding 
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of the interplay between the bulk and deconfinement transitions, we thought 
it necessary to conduct further detailed studies of the £0(3) gauge theory at 
finite temperature. To decouple the two transitions, we 1) employed N T = 2 
lattices in our simulations (a sizeable shift in the transition point was seen 
for N T = 2 lattices in ref. || for the Bhanot-Creutz action) and 2) have used 
extended actions [10] to suppress the bulk transition so as to get a clearer 
view of the deconfinement transition. As in ref. || , we have used the Villain 
form of the £0(3) action. The plan of our paper is as follows: in Sec. |], we 
present results for the phase transition for the N T = 2 and 4 lattices. In Sec. 
[| the phase transition for the extended action is discussed. The last section 
contains a summary of our results and their discussion. 



2 PHASE TRANSITION FOR THE 
VILLAIN ACTION 

2.1 The Action and its Properties 

In the Villain form of the action for the £0(3) theory from ref. B, one in- 
troduces auxiliary Z 2 variables a p associated with the elementary plaquettes, 
besides the usual gauge variables in the fundamental representation of the 
gauge group 577(2). The £0(3) theory is defined by the partition function 

z = E fvUex V (S(U,a)) (5) 
W=±i} 

where the action £ is given by 

S(U,a) = ^Tt f (U p ).tT p . (6) 

On performing the summation over the a p variables, one gets the gauge field 
effective action 

£eff (U) = £ In cosh (yTr/C/^ (7) 

which is clearly blind to the center symmetry. The character expansion of 
e ScS has contributions only from the integer representations of ££7(2) fllQf . 
Explicitly, the link variables have the local Z 2 symmetry 

U t -> -U h a p -> -a p \/p 3 dl. (8) 
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The connection to the SU(2) gauge theory with Wilson action in the weak 
coupling limit is clear: as (3 V — > oo, the configurations that contribute are 

Ut = a ie ^ A , a p = n *i (9) 
ledp 

where o\ are Z2 variables situated on the links. Then a change of variable 
Ui — > Uiai reproduces the SU(2) gauge theory with Wilson action. 

The theory fl5|) is known to have a bulk transition at (3 V ~ 4.45 || |[. 
This bulk transition has been explained by the condensate of Z2 monopoles 
H . The two phases differ in the value of the monopole density, defined by 


M = l-(^E^), (10) 

°c = n °v (n) 

p£<9c 

where c denote the elementary cubes in the lattice and N c is the number of 
such cubes. In the strong coupling phase, 

M — I — 4/3? 

2 

where the character expansion coefficients f3j are given by 

whereas in the weak coupling region, M = up to exponential corrections. 
At the bulk phase transition point, M jumps from its strong coupling value 
to the weak coupling value. 



2.2 Phase Transition at Finite Temperature 

Numerical investigations of the finite temperature transition for this the- 
ory have been carried out in ref. M. By monitoring the plaquette variable, 
P = a p Tif(Up), and (L a ) (Eq. (^)), for lattices with iV T ranging from 4 to 
8, it was found that the theory has only one, strong first order, transition. 
The transition point was found to be f3 v ~ 4.45 for all the lattices, with no 
perceptible change in f3 vc . While this is typical of a bulk transition, the 
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Table 1: Free energy difference for Lee-Kosterlitz analysis. 





6 


8 


10 


12 


16 


24 


AF 


.38(2) 


.52(2) 


.55(2) 


.48(3) 


.58(5) 


.57(10) 



behavior of (L a ) was indicative of a deconfining nature of the transition : it 
jumped across the transition, being ~ before the transition. One possible 
way to explain the curious mixture of the bulk and deconfining characteristics 
of the transition is to hypothesize a shielding of the deconfinement transition 
by a sufficiently strong first order bulk transition. However, since the decon- 
finement transition point is expected to shift with a change in N T , such a 
shielding will only be possible over a finite range of N T ; outside that range 
the two transitions have to separate out. Such a separation was not seen in 
ref. || for N T = 4 to 8. Guided by the results of ref. ||, here we study N T 
= 2 to explore the lower end of the range, as the shift in the deconfinement 
transition point is expected to be larger in going from N T = 4 to 2. 

The transition point was indeed found to shift to (3 V ~ 4.16 for N T = 2. 
The nature of the transition, on the other hand, was found to be different. 
While a two-peak structure was found at the critical point, it did not sharpen 
with an increase in the spatial size of the lattice. We carried numerical 
simulations for lattices with spatial size N a = 6, 8, 10, 12, 16 and 24 (details 
of the simulation are given in Appendix [A]). Figure [l] shows the distribution 
of the plaquette and L a for these lattices at (3 VC . A two-peak structure is 
visible in all of them. However, the peaks move closer to each other as N a 
increases and the valley in between them becomes shallower. This is, of 
course, opposite to what one expects for a first order transition, and suggests 
either a second order transition or a crossover in the thermodynamic limit. 
From Fig. [I] one can calculate the free energy difference between the maxima 
and the minima for a given spatial volume, by using the spectral density 



method J! 1] to extrapolate the histograms to a /3-value such that the two 
peaks are of equal height. For a first order transition, this difference AF is 
expected to increase with lattice size, while a plateau should be reached for a 
second order transition [|l^]. The estimates of AF, shown in table |I|, clearly 



rule out a first order transition. 

Further insight into the nature of the transition can be gained by looking 
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Figure 1: Distribution of a) plaquette and b) L a for N% x 2 lattices with N a 



8, 12, 16, 24. 




Figure 2: a)The plaquette susceptibility, Eq. (fL3|) , as a function of (3 V for N a 
= 6, 8, 10, 12, 16 and 24. The points show the couplings where runs were 
made for each lattice; the lines are Ferrenberg-Swendsen extrapolation, b) 
Peak of the plaquette susceptibility plotted against N c . Also shown is the 
best power law fit, discussed in the text. 
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Table 2: Fit parameters 





a 


b 




x 2 
d.f. 


Xp 


.14 ±.07 


(3.61-?) x 10-2 


1.73 ± .07 


1.12 


Xh a 


-3.6 ±.6 




1.78 ± .07 


.69 



at the plaquette- and L a -susceptibilities 

Xp = 12iV CT 3 «P 2 > - <P> 2 ), (13) 
XL a = N*((Ll) - <L a > 2 ). (14) 

For a first order bulk transition, the plaquette susceptibility is expected to 
scale as the volume of the lattice |13fl . Therefore, for our N T = 2 lattices, 



a 



they are supposed to scale as iV 3 . The peak heights are plotted against N ( 
in Fig. [2] b). Also shown is the best fit to a form XPm.aS^cr) = cl + b N£. 
The corresponding values of the parameters are tabulated in table 0. The 
exponent is clearly different from the expectations for first order transition. 

In case of a deconfinement transition, the peak of the L a susceptibility 
may be expected to scale as iV 3 for a first order transition, and as N£ with 
uo ~ 1.97 for an Ising-like transition. In Fig. |3| b) the peak heights are shown 
as a function of N a . The best fit to the form a + b N£ is also shown. It 
yields u = 1.78 ± 0.07 which, amusingly, is not too far from expectations for 
an Ising-like transition. 

The above behavior is in contrast to what is seen for N T = 4, 6 and 
8 lattices, where a strong first order transition was seen in ref. ||. In 
fact, Fig. [T] is to be contrasted with Figs. 6 and 7 of ref. for N T = 4, 
where the peak positions remained stationary and the peak structure became 
sharper with increase in N a over the range 4 to 8. We have checked that the 
behavior does not alter on increasing the spatial lattice size by making fresh 
simulations for N a = 16 and 24 for iV> = 4. The resulting histograms are 
shown in Fig. [|, along with the 8 3 x 4 data from ref. [Q. They clearly 
indicate a first order transition for N T = 4 lattices in the infinite volume 
limit, with discontinuities in plaquette and adjoint Polyakov loop given by, 
AP = .058 ± 3, AL a = 0.87 ± 0.04. 

The second order transition for the N T = 2 case, with an exponent close 
to the 3-D Ising model, may prompt one to identify the transition with the 



7 





6 — 


j 


8 
10 


- a) 1 


12 


i 16--- ' 




24 - - 

% 




Hi 

— \ Nqt^" I. _ 








Figure 3: a)The Polyakov loop susceptibility, Eq.([I4|), as a function of (3 V for 
= 6, 8, 10, 12, 16 and 24. b) Peak of the Polyakov loop susceptibility 
plotted against N a . Also shown is the best fit discussed in text. 




Figure 4: Distribution of a) plaquette and b) L a for N% x 4 lattices with N a 
= 8, 16, 24. 
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SU(2) deconfinement transition, thereby vindicating the shielding scenario 
above. However, on a closer examination, one sees that the detailed proper- 
ties of the transition are very different from the SU (2) case. The plaquette 
variable is smooth at the SU(2) deconfinement transition point, and \L a does 
not show criticality there ||14|| , whereas from table |^ one finds that both \L a 
and xp diverge at the SO (3) transition with essentially the same exponent. 
This indicates that here too one is seeing an intricate interplay of bulk and 
deconfinement transition effects, very much like the higher N T cases. We be- 
lieve that the reason the order of the transition changes in going from higher 
N T to iV T = 2 may simply be the fact that the latter have effectively a smaller 
bulk-dimensionality. This interpretation, suggested by the near equality of 
the exponents for \L a and xpi indicates a more complicated interplay of the 
bulk and deconfinement transitions than the simple shielding scenario since a 
shielding is easy to justify only in the vicinity of a strong first order transition, 
while the N T = 2 lattices show a mixing of effects of bulk and deconfinement 
transitions even in the presence of a continuous bulk transition. Moreover, 
it also calls for a different approach to disentangle the bulk effects from the 
deconfining ones. In the next section, we present results from an attempt in 
that direction. 

3 PHASE TRANSITION FOR 
EXTENDED ACTION 

Since the SO (3) bulk transition is known to be caused by a condensation 
of Z 2 monopoles, which are purely topological lattice artifacts and do not 
survive in the continuum limit, one may hope to get a clearer view of the 
physics of the deconfinement transition in the SO (3) theory by suppressing 
them, without missing any important continuum physics. One may be able 
to address important physics questions associated with the deconfinement 
transition in this manner, specially since there is no obvious order parameter 
in this case. The suppression can be achieved smoothly by adding a term to 
the action which disfavors the monopole configurations. The most obvious 
choice is to add a monopole chemical potential term to the action [|H]], as 
considered below. 
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3.1 Action and its properties 



We study the action [JlQ~ | 



Sm{U, cr) 



Pv 



(15) 



where a c is defined in Eq. (|TT|). 




Figure 5: The phase diagram for action flip]). The dotted line denotes a sec- 
ond order transition while the dot-dashed line denotes a first order line. The 
points are the actual simulation points. The full line is the approximation of 
Eq. (H). 



This action has a line of phase transitions in the (/3 V , A ) plane [ 10| . For 
A = 0, it reduces to the action of Eq. ([]), with a first order transition at 
j3 vc ~ 4.5. For (3 V = 0, the theory is dual to four dimensional Ising model 
with coupling A = |ln coth A |T5]] . It therefore has a second order phase 
transition, at A c ~ .95. For small (3 V , one can integrate out the gauge field 
: this contributes, in the leading order, to only a renormalization of the 
coupling 

A — > \ eff = A + 4/3| (16) 

where (3i is given in eq. ([12]). Therefore, the transition will remain second 
order, at least to a point where the leading order approximation is valid. As 
j3 v is increased, it becomes first order and joins the transition point on the 
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(3 V - axis. We show the bulk phase diagram in Fig. Apart from some 
details like the second order line etc., it is qualitatively the same as that in 
ref. |[L0|| . For large /3 V , the large A region is smoothly connected to the small 
A region. For, due to Eq. (J^) the second term in action fll5|) is irrelevant 
here (any other configuration is strongly suppressed and can be treated in 
the dilute gas approximation). So the weak coupling limit is not changed by 
the addition of this term. 

3.2 Study of the Finite Temperature Phase Transition 

In order to study the true deconfinement phase transition which is unaffected 
by the bulk effects, we need to choose a value of A > A c , so that the bulk 
transition can be avoided. For our simulations, we have used A = 1. We 
checked that with this value of A the plaquette variable is indeed smooth ev- 
erywhere, indicating the lack of any bulk phase transition. The approximate 
location of the deconfinement transition is indicated by a steep rise in (L a ) . 
We studied the deconfinement transition for lattices with N T = 2 and 4. For 
finite size analysis, lattices with N a = 8, 12 and 16 were used for both N T . 

For N T = 2, the deconfinement transition was found to be at (3 V ~ 1.9, 
while it was found to be at (3 V around 2.3 for the N T = 4 lattices. Both of 
these are close to the corresponding transition points for pure SU{2) and 
far from the transition points for action (^j) (see Sec. |2.2|) . We carried out 
simulations at a few f3 v around the approximate transition point for each 
lattice size and then used multi-histogramming methods |1(J to obtain the 
physical variables in the whole critical region. The details of the simulation 
are given in Appendix In Fig. |6] the behavior of (L a ) across the transition 
point, as obtained by the multi-histogramming extrapolation, is shown. (L a ) 
is seen to rise smoothly without any signs of a discontinuity in it. This is 
similar to its known behavior across the SU(2) deconfinement transition. 
We also found that the the (L a ) susceptibility did not show any peak at the 
transition point, again in agreement with the expectations from the SU(2) 
theory. 

A quantitative study of the order of the deconfinement transition is made 
difficult by the absence of an order parameter for this transition. The two 
physical variables that can be used are the energy density and specific heat. 
Extraction of these observables from lattice simulations is, however, difficult, 
as they are not simply related to objects directly measured in simulations. 
Since our purpose is only to show the existence of a deconfinement transition, 
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Figure 6: a) (L a ) across the deconfining point for N T = 2 lattices with N a 
= 8, 12, 16. The points with error bars show the actual couplings where 
simulations were made. The curve is obtained by using multi-histogramming 
methods, b) Same for N T = 4. 



and not to make quantitative predictions about the transition, we use for 
these variables the approximate expressions (see Appendix |B]) 

ea 4 = 3p v (P T -P ), (17) 
O 3 = 3N T N p (3 2 var(P T -P a ), (18) 

where N p is the number of spatial plaquettes, P T and P a are spatial and 
temporal plaquette variables respectively, and var(x) = (x 2 ) — (x) 2 . The 
energy density for N T = 2 and 4 lattices is shown in Fig. [F[ It is seen to 
rise significantly across the transition for both the N T values, justifying the 
identification of the transition as a deconfining one. The rise becomes steeper 
with increase in the spatial volume, indicating a true transition as opposed 
to a crossover. 

The behavior of the specific heat for the N T = 2 lattices is shown in Fig. || 
The position of the peak corresponds to the critical point. The peak height 
is seen to increase with the lattice size, and the peaks become sharper. Both 
these behaviors are expected for a deconfinement phase transition. While 
our expressions for specific heat, Eq. fll~8|), are somewhat approximate to 
obtain an accurate critical exponent, it is interesting to note that a linear 
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Figure 7: The energy density, Eq. (|I7D , in the critical region of a) N% x 2 
and b) iV ff 3 x 4 lattices for N a = 8, 12, 16. 



fit to the peak heights gives ^ = 0.13 ± 0.04, which is consistent with the 
three dimensional Ising model value - ~ .17. For the N T = 4 lattices, the 
energy density being 16 times smaller in lattice units, the resultant specific 
heat curves were very noisy and inconclusive; a measurement of the exponent 
here is computationally very demanding and beyond our reach. 

Thus the deconfinement transition for the monopole suppressed action 
of SO (3) is very similar in nature to that for the SU(2) Wilson action, 
although the former has no order parameter. A significant difference in the 
two, however, is the presence of negative (L a ) states in the SO (3) theory. 
For the SU(2) Wilson action, (L a ) approaches a unique positive value — > 
3 as (3f — > oo in the high temperature phase. However, for 5*0(3), (L a ) 
manifests itself in two states in the high temperature phase. One is like the 
high temperature state of SU(2) but the other has a negative (L a ) , which 
approaches —1 as (3 V — > oo. In ref. ||, it was proposed that the negative (L a ) 
state corresponds to a new bulk phase which is separated from the (L a ) ~ 
bulk phase by a bulk phase transition. From the behavior of the correlation 
functions, and the behavior of (L a ) in presence of a "magnetic" term hL a in 
the action, it was concluded in ref. || that both the positive and negative 
(L a ) states correspond to the same physics. In the monopole-suppressed 
action, there is no bulk transition. However, we still find the negative (L a ) 
state, very similarly to the monopole unsuppressed action. In Fig. [5] we show 
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the behavior of (L a ) for a large range of the coupling (3 V , for both A = 
(the pure Villain action) and A = 1 (the monopole suppressed case). Other 
than the shift in transition point, there is no essential difference between the 
two. Therefore the possibility advocated in ref. || of the negative (L a ) state 
being another bulk phase, separated by a bulk transition, is disfavored. 

4 Summary and Discussion 

We studied the 5*0(3) lattice gauge theory in a Villain form numerically. Our 
main objective was to disentangle the effects of the bulk phase transition in 
this theory from those of the deconfmement phase transition in it, if any, and 
then to compare the latter with the SU(2) theory which is expected to be in 
the same universality class. We attempted this by investigating the case of 
smaller N T = 2 lattices and later by suppressing the Z 2 monopoles. Going 
to N T = 2, the nature of the only transition observed was seen to change 
substantially compared to the earlier N T = 4-8 studies. The transition point 
had a visibly large shift, and the order of the transition was found to change to 
second. However, the curious mix of bulk and finite temperature transitions, 
seen on lattices with higher N T , continued, as the behavior of the plaquette 
variables and that of (L a ) , the adjoint Polyakov loop, were very similar. 
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Figure 9: (L a ) for a range of couplings, showing the two different states in 
the high temperature phase. The filled circles and triangles correspond to 
the Villain action of Sec. |2] while the open circles and triangles correspond 
to the monopole suppressed action. 



Determination of the critical exponents for their respective susceptibilities 
suggested a coincident bulk and deconfinement phase transition, similar to 
the higher N T studies but with a second order bulk transition. The scenario of 
bulk transition shielding the deconfinement transition, which was a plausible 
one for the higher N T lattices, is therefore disfavored, indicating the necessity 
of a more complicated scenario. 

In order to get a clearer view of deconfinement in the 5*0(3) theory, the 
bulk transition was evaded by adding a term in the action to suppress the 
topological objects driving the bulk transition. The resulting theory was 
found to have a deconfinement transition whose nature is very similar to the 
deconfinement transition in SU(2). This is very interesting as the center 
symmetry, which plays a crucial role in the SU(2) deconfinement transition, 
is trivial here, and the theory does not have any order parameter for its 
deconfinement transition. The location of the transition is also seen to be 
very different from that where the coincident bulk/deconfinement transition 
was seen for the Villain action without monopole suppression. It seems, 
therefore, that the bulk transition affects the deconfinement transition in 
a very subtle and nontrivial way, making it always coincide with itself. It 
would be very interesting to see two separate transitions for the original 
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SO (3) theory without any monopole suppression on a big enough lattice, 
since these monopoles are mere lattice artifacts. 



5 Acknowledgments 

One of us (S.D.) would like to thank Profs. Frithjof Karsch, Sumit Das and 
Sourendu Gupta, and the other (R.V.G.) would like to thank Dr. Manu 
Mathur for stimulating discussions and suggestions. 



A Details of Numerical Simulations 

In this appendix we give details of the numerical simulations mentioned in 
the text. For all the simulations, we have used heat bath algorithms for 
both the gauge variables and the a p . An iteration consisted of updating all 
the link variables once, followed by an update of all a v . For the N T = 2 



measurements of Sec. 2.2, the transition point was approximately located 



for each lattice by looking for a phase coexistence. At this point a small run 



was made, and the spectral density method |TT| was used to locate the sus- 
ceptibility peak where a long simulation run was made. We used 1.6 million 
configurations for each lattice. Table |] lists the couplings at which the long 
runs were made for different N a and the autocorrelation times for plaquette 
and L a . For the extrapolations to nearby couplings using the spectral density 
method, the Jackknife estimate was used for both the extrapolated value and 
the error, using 20 blocks for the analysis. For the histograms of 16 3 x 4 and 
24 3 x 4 lattices of this section, we produced 50000 configurations for each 
lattice, which was sufficient to clearly establish the peaks and nature of the 
histograms. The phase coexistence region was much broader for these lat- 
tices; we took (3 V = 4.45, approximately at the middle of the region, for the 
longer runs. Due to the strong first order nature of the transition, tunnelings 
are not seen and we had to combine configurations from two runs to get the 
peaks. 



For the simulations in Sec. |3]2|, it was noticed that the plaquette variables 
have an extraordinarily large autocorrelation time. This can be reduced dra- 
matically by implementing the energy- conserving steps corresponding to the 
center symmetry transformation (Eq. (|8])). After every combined step of the 
previous section, we added another sweep where a fixed fraction (arbitrarily 
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Table 3: Details of the Monte Carlo simulations of Sec. 



2.2 





Pv 


rp 


TL a 


6 


4.16 


516.8 


754.7 


8 


4.159 


1010.4 


1294.3 


10 


4.156 


1738.5 


1878.8 


12 


4.156 


2245.7 


2254.6 


16 


4.1552 


2477.6 


3896.5 


24 


4.1548 


3155.8 


4753.1 



chosen to be one-fourth) of the total number of links, randomly chosen, are 
flipped along with all the plaquettes touching the link. An iteration con- 
sisted of the combination of this sweep with the combined heat bath step, 
and measurement was taken after every iteration. It was checked that for all 
the runs, the autocorrelation time was less than 100. For the N T = 2 results 
of this section, from the rise in (L a ) the critical point was estimated to be 
~ 1.90. We ran simulations at (3 V = 1.82, 1.85, 1.87, 1.88, 1.90 and 1.92, 
producing half million configurations at each coupling for each lattice. For 
the N T = 4 case, the critical point was at ~ 2.30. We produced one million 
configuration at j3 v = 2.30, and half million configurations each at (3 V = 2.25, 
2.28, 2.32 and 2.35. For interpolating to other values of (3 V , each dataset 
was divided into very finely meshed histograms, and the multihistogramming 
techniques of ref. |16| was applied over the whole set. The central values 
and the errors quoted are Jackknife estimates, 20 blocks being used for the 
analysis in each case. 



B The Energy Density and Specific Heat 

The expressions for energy density and specific heat density, Eqs. ( |17D and 
(|j~8|), can be easily obtained using the methods of ref. [^7f| . We outline here 
the derivation of these results, emphasizing the approximations involved. 

For finite temperature theory, it is necessary to take into account the 
asymmetry in the spatial and temporal directions by using two different 
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couplings. Writing 

S = p a £ P„ + (3 T £ P T , (19) 

spatial p temporal p 

the correct continuum limit is recovered by writing 3 a = 4-— and /? T = 4- — , 
where a T and a CT are the temporal and spatial lattice spacings, respectively. 
The couplings g a and g T both equal the usual bare coupling g when the two 
lattice spacings are equal. In general, one can expand them perturbatively 
in g near the continuum limit [|T7|]: 

g~ 2 {a a , a T ) = g~ 2 (a a ) + c ajT (a a /a T ) + 0(g 2 ). 

Defining the inverse temperature (3 = N T a, the energy density is given by 

1 d 



vdd AnZl 



N^alN T \da T la 



(20) 



From Eq. ([19]), one sees that evaluation of dS/da T \ aa requires knowledge 
of the nonperturbative beta functions dg~^/da T \ atT , which, however, are not 
known for 5*0(3) gauge theory. One can, therefore, use the perturbative 
expressions for the same, using the expansion in terms of g~ 2 given above. 
The perturbative beta functions are exactly similar to the SU (2) case of ref. 



171 . The contribution from this term is found to be small and does not 
affect the critical behavior. Ignoring this term, and putting 
after taking the derivatives in Eq. (p0| , one gets the expression (|T7|) for the 
energy density. 

Similarly the specific heat at fixed volume is 

= ^{2f3N p {P T -P )+f3 2 N 2 p ^i{P T -P a )} 1 (21) 

where P D is the plaquette variable on a symmetric lattice, and the coupling 
derivative terms have again been neglected. The first term in Eq. ( pi] ) 
is of same order as the neglected derivative terms, and should be omitted 
for consistency, giving us the expression (|T8|). We have also checked that 
inclusion of the first term does not change the results in any significant way. 
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